Entanglement of electron spins of non-interacting electron gases 
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We study entanglement of electron spins in many-body systems based on the Green's function 
approach. As an application we obtain the two-particle density matrix of a non-interacting electron 
gas and identify its two-spin density matrix as a Werner state. We calculate entanglement measures, 
a classical correlation, mutual information, and a pair distribution function of two electrons at zero 
and finite temperatures. We find that changes of entanglement measures are proportional to T 2 at 
low temperatures. 



PACS numbers: 03.67.Mn, 03.67.-a, 03.65.Ud, 71.10.Ca 



Introduction.- Entanglement is considered to be one 
of the key resources in quantum information science P, 
0, . Much attention has been paid to quantifying de- 
grees of entanglement, to generation of entangled states, 
and to applications of entangled states to quantum com- 
munication and quantum teleportation |2J, y| . Recently, 
considerable interest has been devoted to entanglement 
of two subsystems of a many-body system: quantum 
spin wstemsjljljl 0, H 03, identical parti- 
cles FLU [H Ellillia 0, fractional quantum Hall ef- 
fect ^3 1 an d spins of a non-interacting electron gas . 
Entanglement shows non-classical (or nonlocal) correla- 
tions between quantum systems. In many-body systems 
the correlation functions play a fundamental role in de- 
scribing their physical phenomena. Thus it is natural to 
explore the relation between entanglement and the cor- 
relation functions. 

Since up to now entanglement measures have been rel- 
atively well developed for two qubits, one needs a two- 
particle density matrix to study entanglement of a many- 
body system. If a state of a many-body system is known, 
then it is possible to calculate its useful physical quan- 
tities and also a two-particle density matrix directly by 
tracing out the rest of the system. However, it is impos- 
sible to obtain exact many-body states except for very 
simple systems. Instead of finding a many-body state, 
usually one works with Green's functions which make it 
possible to study the effects of interaction in a system- 
atic way. In this paper we adopt this way to study the 
entanglement of many-body systems. As an illustration, 
we investigate the entanglement of two electron-spins of 
free electron gases at zero and finite temperature follow- 
ing Vedral's work 0. We find the two-spin state of a 
non-interacting electron gas is given as a Werner state. 
We also discuss the relation between entanglement mea- 
sures, classical correlations, the total correlation, and the 
pair distribution functions. 

Entanglement of two- electron spins.- Consider a sys- 
tem of N non-interacting electrons in a box with volume 
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V. The ground state of the system is 

i*o>= n c Lio), 

|k|<fc F 



(1) 



where kp — (S^N/V) 1 ^ 3 is the Fermi momentum. From 
the density matrix = |*o)(^ r o| of the system, it is 
easy to obtain the two-particle density matrix |23| 



p {2) {xi,x 2 ;x' 1 ,x' 2 ) = ^ 



p^{x 1 ;x' 1 ) p^{x 1 ;x' 2 ) 
pW(> 2 ;xi) p^(x2;x' 2 ) 



(2) 



where p {1 \x;x') = J2k 4>t<j( Y )4>\t<j' ( r ') is the one- 
particle density matrix. Here x = (r, a) denotes the posi- 
tion and spin quantum numbers of an electron, <f>^ a (r) = 



e Xo-j an d Xcr is the spin wave function. 



The Green's function approach is very convenient in 
solving many-body problems. The two-particle density 
matrix is given by 

p^(x 1 ,x 2 ;x' 1 x' 2 ) = , (3) 

where (O) = (^ol^l^o) f° r zero temperature and (O) = 
Tr{pcO} for finite temperatures, T / 0, with Zq = 



Tr{e-ft H -i* N )} and p G = e 



/Zq. The two- 



particle temperature Green's function is defined by 
0(1, 2; 1', 2') = Tr{p G T T [^(l)fe(2)^(2')^(l')]},(4) 
where the number 1 denotes the variable (xi,t±). The 



field operator is defined by ipx (xt) 



^Kt/h 



ip{x) e 



-Kt/H 



with K = H — fiN. The relation between p^ and 
0(1,2; T,2') is given by 
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p { >(x 1 ,x 2 ;x' 1 ,x 2 ) = --g{x 1 T 1 ,x 2 T 2 ;x' 1 T' 1 + ,x' 2 T 2 + ) , (5) 

where t + denotes a time infinitesimally later than r. 

In general, it is difficult to find the exact Green's func- 
tion for an interacting many-body system. One of the ap- 
proximations to 0(1, 2; 1', 2') is the Hartree-Fock approx- 
imation. A generalized Wick's theorem makes it possible 



2 



to express the two-particle temperature Green's function 
in terms of one-particle Green's functions approximately 

0(1, 2; 1', 2') a 0(1, l')G(2, 2') - 0(1, 200(2, 10 • (6) 

where 0(1; 10 = Tr{p G T T [^ K (l)^(lO]} is the one- 
particle temperature Green's function. Beyond the 
Hartree-Fock approximation of Eq. ©, one needs the 
calculation of the vertex part [2l| • 

For a non-interacting system considered here, Eq. © 
is exact. Also, it is easy to calculate the non-interacting 
Green's function 0°(1; 10 from its definition or by con- 
structing an equation of motion for 0°(1; l0< One obtains 



pW(x; x 1 ) = -G°(xr; x't+) = <W g(r - r') , (7a) 
where the one-particle space density matrix g(r—r') reads 
1 



7(r) 



k 



(7b) 



Here rik = {exp[/3(ek — + is the mean occupation 
number in state k with energy = h 2 k 2 /2m. At zero 
temperature one has = 9{kp — |k|). 

With Eqs. (JjjJ, JJjJl, and J2J), one has the explicit form 
of the two-particle space-spin density matrix |19l |23j 

i( 2 ) 



p ( ) {xi,x 2 ;x' 1 ,x l 2 ) 

= g[ f( r i _ r 'i)ff( r 2 - r 2 )<Wi<W^ 
- -r 2 )ff(r 2 -ri)<5 CTlCT /5 CT / CT2 ] 



(8) 



To the best of our knowledge, it seems that there is no 
entanglement measure of identical particles, which takes 
into account both continuous variables and discrete in- 
ternal variables. Depending on the space density ma- 
trix, two spins may be entangled. We obtain Vedral's 
result 18] only if ri = and r 2 = r 2 , that is, only diag- 
onal elements of a space density matrix are considered. 
The two-spin density matrix, depending on the relative 
distance between two electrons r = \vi — r 2 |, reads 

P ( 'L-y^(r) = f [<W(<W 2 - /(r) 2 <W'^ CT2 ] , (9) 



where n = N/V is the particle density and 
/(ri - * J- - ■ \ » 



(10) 



In condensed matter physics, correlation between spins is 
described by two pair distribution functions, = gn = 
(1 - / 2 )/2 and g n = 1/2 ^ The integration relation 
J f 2 dr 3 = 2V/N gives rise to the normalization condi- 
tion of p {2 \ Tr-^ 2 '} = N(N-l)/2. At zero temperature 
the analytic form of f(r) depending on the spatial dimen- 
sion of the system is well known 

3ji(kpr)/kpr , 3-dimension 

m = { , (ii) 

2 J\(kpr) /kpr , 2-dimension 



where j\ is the spherical Bessel function and J\ the 
first-order Bessel function of the first kind. Here 2- 
dimensional Fermi wave vector is kp = \f2wn. 

By dividing the bracket part of Eq. © by 4 — 2/ 2 , we 
get the two spin-density matrix p\ 2 for a given relative 
distance r between two electrons 



Pl2 = 





rw 2 








o ■ 


1 





1 


-p 





4-2/2 





-f 2 


1 
















1-/ 2 J 



(12) 



where Tr CT1<T2 { j o 12 } = 1. We find that p 12 is nothing but a 
Werner state characterized by a single parameter p [24| 

I 



P12 = (1-P) 7 +P|*~>(* !■ 



(13) 



where I is a 4x4 identity matrix, |* _ ) = (|01)-|10))/\/2, 
p = f 2 /(2 - f) with < p < 1, and the fidelity F = 
(*-|/> ia |*-) = (3p + l)/4 = (/ 2 + l)/(4 - 2/ 2 ) with 
1/4 < F < 1. It should be noted that a Werner state 
appears in an anti- ferromagnetic Heisenberg model Q. 
In fact the exchange interaction between free electrons 
gives rise to the anti-ferromagnetic coupling. 

The properties of a Werner state are well studied. Ac- 
cording to the separabilit y cr iterion of two qubits by the 
partial transposition H^Jl^, if p > 1/3 (i.e., F > 1/2 
or f 2 > 1/2), then p\ 2 is entangled. For p > 1/V2 (i.e., 
F > (2 + 3\/2)/8 or f 2 = 2( v / 2 - 1)), p u violates the 
Bell-CHSH inequality |2jfl. There are a few computable 
entanglement measures for two qubits. The concurrence 
for p 12 [H is given by C = max{0, (2/ 2 - l)/(2 - f 2 )}. 
Using C, one has the entanglement of formation Ep(C) — 

h( 1+v/ ^~ — ) where h(x) is the Shannon entropy. For a 
Werner state, one has the relative entropy of entangle- 
ment E RE {pi2) = 1 + Flog 2 F + (1 - F)log 2 (l - F) 
for 1/2 < F < 1 IH Hg. We take E RE (p 12 ) = for 
1/4 < F < 1/2. For 1/2 < F < 1, is expressed in 
terms of / 

i + / 2 1 , r i + / 2 



E RE (pi 2 ) = 1 + 



+ 3 



4-2/2 



log 2 



2/2 



log 2 



-2/2 

I-/ 2 



2/2 



.(14) 



Recently, some attention has been paid to the split- 
ting of classical and quantum correlations from the total 
correlation J^l], [3^, |33 ■ Usually the total correlation is 
given by the mutual information 



I{pi2) = S(p!) + S{p 2 ) - S{p 12 ) 



(15) 



where S(p) — — Tr(plog 2 p) is the von Neumann entropy. 
For the Werner state one easily obtains I(p\ 2 ) = 2 + 
Flog 2 F + (1 - F) log 2 [(l - F)/3], which can be written 
in terms of / 



+ 3 



1 + / 2 
4-2/2 

I-/ 2 
4-2/2 



log 2 
log 2 



1 + / 2 

4-2/2 



4 - 2/2 



(16) 
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It should be noted that Eq. JTJJ differs from Eq. (47) 
in Ref. 0. If f 2 = 0, 1 then I(p 12 ) = 0, 2, respectively. 
Among definitions of the classical correlation |3ll l32l l33| . 
we follow Hamieh et al.'s definition of the classical cor- 
relation by C c /(pi2) = I(pi2) — Ere 33]. For the 
Werner state, we have C c i(pi2) = 1 — (1 — F) log 2 3 for 
1/2 < F < 1 and C c i(p 12 ) = I(pi 2 ) for 1/4 < F < 1/2. 
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FIG. 1: The mutual information I(pi2) (thick solid line), the 
relative entropy of entanglement Ere (dashed line), the clas- 
sical correlation C c ; (dotted line), the concurrence C (solid 
line), and the function f 2 (dashed-dotted line) as functions of 
fcyr for a 3-dimensional electron gas at zero temperature. 

Fig.^shows the relative entropy of entanglement Ere, 
the concurrence C, the mutual information /(/O12), the 
classical correlation C c ;(pi2), and function f 2 as a func- 
tion of the relative distance kpr of two electrons normal- 
ized by the Fermi wave length 1/kp at zero temperature 
for a 3-dimensional electron gas. The shapes of all the 
functions for a two-dimensional electron gas are similar 
to those for a 3-dimensional gas. Usually it has been 
known that the more noticeable oscillation of /(r), the 
stronger the exchange correlations of electrons. However, 
the entanglement measures, Ere and C, show no oscilla- 
tory behavior. We see that the behavior of the classical 
correlation C c i is similar to that of f 2 . Unlike the con- 
currence C, the relative entropy of entanglement Ere is 
differentiable at f = 1/2 (F = 1/2). 

As shown by Vedral Jig, one can expect the entan- 
glement of two spins within the order of the Fermi wave 
length l/kp. In usual metals, the Fermi wave length is 
the order of A. However, we would like to point out that 
the Fermi wave length of a 2-dimensional electron gas is 
the order of hundred A. Thus it may be possible to ex- 
tract entangled spins out of a 2-dimcnsional electron gas 
formed in GaAs heterostructure. 

In order to study the entanglement at finite tempera- 
ture, we should evaluate Eq. (|10fl rewritten by 



f{r,T) 



2k F r 



s\a.(kpr^/x )n{e F x)dx , (17) 



where n(ek) = and the argument T is explicitly shown 
in order to emphasize the temperature dependence of /. 



We calculate Eq. I|17fl by the numerical integration and 
by the Sommerfeld expansion. Fig. [5] shows f(r,T) and 
Ere at two normalized temperatures, T/Tp = and 
0.15, where Tp = ef/^b is the Fermi temperature. Fig. [3] 
plots A/(r,T) ee f(r,T) - /(r,0) where /(r,0) is given 
by Eq. l|TT]l. 
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FIG. 2: (color online) Relative entropy of entanglement Ere 
and / as functions of kpr at temperatures T/Te = (solid 
line) and 0.2 (dashed line) for the 3-dimensional electron gas. 
The inset shows the mean occupation number at corre- 
sponding temperatures. 
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FIG. 3: (color online) A/(r, T) as a function of k F r and T/T F - 



Using the Sommerfeld expansion, we see the temper- 
ature dependence of f(r,T) in more detail. At kpr ~ 1 
we obtain the approximation of A/(r, T) with the order 
of T 2 as 



A/(r,T) 



^2 r 



cos(kpr) 



sm(kpr) 



kpr 



k B T 



(18) 



Fig- El ( a ) shows the plots of A/(r, T) calculated by the 
numerical integration of Eq. (17) and the right hand 
side of Eq. i|18[l obtained by the the Sommerfeld ex- 
pansion. We see that at low temperatures the shift 
of / is the order of T 2 . The change of f(r,T) gives 
the change of entanglement measures. We calculate 
AC(r,T) ee C(r,T) — C(r, 0) written approximately as 



AC(r,T) 



6/(r,0) 



[2-/(r,0)2 



fA/(r,T) 



(19) 
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FIG. 4: (color online) (a) Af(r, T)/(T/T F ), (b) 
AC(r,T)/(T/T F ), and (c) AE RE /(T/T F ) versus T/T F 
at fe_FJ" ~ 1. The dashed lines are obtained from the Som- 
merfeld expansion. The solid lines are due to the numerical 
integration of Eq. 1171 . 



where in derivation we keep only the terms of order 
A/(r, T), Thus AC is proportional to T 2 at low tem- 
peratures as depicted in Fig. 0] (b). Also we evaluate 
AE RE (r,T) = E RE {r,T) - E RE {r,0) and find the T 2 
dependence of AE RE at low temperatures as shown in 
Fig. E|(c). 

Summary.- In this paper, we obtained the two- 
particle density matrix of a non-interacting electron gas 
based on the Green's function method. It was shown 
that the two-spin density matrix for a given relative dis- 
tance between two electrons has the form of a Werner 
state of which the parameter p is a function of the space 
density matrix /. We presented the relation between the 
total correlation, the entanglement measures, the classi- 
cal correlation, and the pair distribution functions. Also 
we have shown that the entanglement measures change 
a little bit in proportion to T 2 at low temperatures and 
fcj?r ~ 1. 

Some remarks should be made. First, in discussion of 
entanglement of two spins we ignored the space density 
matrix of electrons. It is interesting to find entanglement 
measures for identical particles with spatial and internal 
degrees of freedom. Second, in this work the electron- 
electron interaction was ignored. It will be our future 
work to investigate how the electron-electron interaction 
influences the entanglement measures of electron spins. 
Finally, how to extract entangled spins out of an electron 
gas would be another interesting problem. 
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